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Bell polynomials basis. The properties of Bell polynomials are explained. We 
also introduce the operational matrix of derivative for Bell polynomials. The 
chief feature of this matrix is reducing the OCPs to an optimization prob- 
lem. Finally, we discuss the convergence of the new technique and present 
some illustrative examples to show the effectiveness and applicability of the 
proposed scheme. Comparison of the proposed method with other previous 
methods shows that this method is accurate. 
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1 Introduction 


Researchers believe that the analytical solution of optimal control problems 
(OCPs) is usually complicated, or, most of the time, impossible [14]. Hence 
we use numerical methods to solve them. The first attempt to solve OCPs was 
applying numerical methods to related Bellman and Ponteryagin in 1950; see 
[16]. The main problem in numerical methods is the accuracy of computations 
and the amount of manipulations. 


In recent years, some studies have been done in this regard, which we 
will refer only [1, 7, 8, 9, 11, 18, 20, 21, 25]. The investigation of OCPs and 
solving them by using the Bezier polynomials were proposed by Yari, Mirnia, 
and Lakestani [25]. To be more specific, they chose Bezier polynomials as 
basic functions and considered the state vector «(t) and control vector u(t) 
as Bezier polynomials with unknown coefficients. Grigoryev, Mustafina, and 
Larin [9] proposed the numerical solution of OCPs, using a successive approx- 
imations method. Ramezani [18] presented the numerical solution of OCPs, 
using the second kind Chebyshev wavelet. Rose [20] described numerical 
methods to resolve control issues. Kafash and his colleagues [11] proposed a 
numerical approach for solving OCPs, using the Boubaker polynomials ex- 
pansion scheme. 


In the present study, we offer a numerical solution for the linear con- 
strained quadratic OCPs via Bell polynomials basis. The Bell polynomials 
were studied extensively by Bell [4] in 1934. These polynomials naturally oc- 
cur from differentiating a composite function several times. Bell polynomials 
have many applications in number theory and classical analysis, and there 
is a vast literature about their applications. They are frequently applied in 
combinatorial analysis and statistics. Also, these polynomials have been used 
in many other contexts such as the Blissard problem; see [15]. 

This method relies on approximating the state and control variables with Bell 
polynomials. It includes reducing the OCPs to an optimization one by first 
expanding the state rate «(t) and the control u(t) as a Bell polynomial with 
unknown coefficients. The operational matrix of differentiation Dg is given to 
approximate the differential part of the problem and the performance index. 


This paper is arranged as follows. In Section 2, we present some important 
preliminaries needed for our subsequent development and approximation of 
a function by using the Bell polynomials basis. Section 3 is devoted to the 
numerical method of solving the OCPs. In Section 4, we discuss the proposed 
method for solving the OCPs. In Section 5, we present the convergence of 
the proposed method. In Section 6, we consider the applicability as well as 
the accuracy of the proposed technique with several examples. Finally, the 
main results are summarized in Section 7. 
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2 Preliminaries 
In this section, we present some preliminaries needed in some other parts of 


the paper. 


2.1 Bell polynomials and their properties 


The incomplete exponential Bell polynomials are given as 


m 1 dy ( tm—n+1 )dm—n+1 , 


Brn (#1) -+-1tm—m+1) = lade! i GaSe PI) 


m™(m,n) 


where m,n € Zt, n < ‘mz, and the summation takes place over all integers 
dy, da, ebay dm—n+1 = 0, such that 


dy + dg+-++++dm-n4i =, 
dy + 2d2 +---+(m+n—1)dm—en4i =m. 


The sum se 


Brlbu trast => Beal tte ta ea); 


n=1 
is called the mth complete exponential Bell polynomial [4]. 
Suppose that t; = tg = --- = tm = t, and let B,,(t,t,...,t) = Bn(b), 
then 


Bm(t) = So {m,n}t”, (1) 
n=0 
is called the complete Bell polynomial (Touchard polynomial), where 
Le 3 (nr age 
{m,n} = — > (-19(" \n—)", 
ie j 
j=0 
are called Sterling numbers of the second kind and some properties of them 
are as follows (see [4, 23]): 
1) {m,0} = bm, 
2) {m, 1} =1, 
3) 4,2) 2"! 1, 


4) {m,n} = n{m—1,n}+{m-—1,n—-]}, 


5) {m,m} = 1, 
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(6) {m,n} =0,n>m. 


Some sterling numbers of the second kind are given in Table 1: 


Table 1: Sterling numbers of the second kind 


{m,0} {m, 1} {m, 2} {m, 3} {m, 4} {m,5} {m, 6} {m, 7} {m, 8} {m, 9} 
1 
3 1 
4 6 al 


15 25 10 1 

31 90 65 15 1 

63 301 = 350 140 21 1 

127. =966 =: 1701 )=1050) =. 266 28 1 

255 3025 7770 6951 2646 462 36 il 


WOON Ao RwNMr Oo] 3 
eo cS © 2 oo 2 oO © = 
Pe PRP RP RP eee 


24 


(t) 

(t)=t 

(t)=t 
Bs(t) = ¢° + 3¢? +4, 

(t)=t 

(t)=t 

(t)=t 


Phe te Ask 
Ph Oe Dae 15e? ee, 
6 4 15 + 65 + 90P + a1 +e 


Some properties of Bell polynomials are as follows (see [4, 5, 6]): 
(1) By(1) = By (Bell number), 
(2) Bm(—1) = Bm (complement of Bell number), 
(3) for all m,n € Z* and n<m, 


I, 


Je at -oaeg Jie at _ B;(at)e”, 
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: Bini (t) = t(Bm(t) + Bi, ()), (2) 
4. d 
tT Bm(te! = Bmsi(tye’, 
5. 3 
t 7 Bm( te" =t(Bm(t) + Bi (@))e’, 

6. 2 

Bmarlt)=t)0 (™) Balt (3 

n=0 


m 


Birla) => (") Ba(P)Bmn-n(0) 


n=0 
Theorem 1.|semi-orthogonality of B,,(t)] For all m,n € N, we have 
e72t gm+rn a | 


= yyr-t 
, ( ) mtn 


i Bm ( BA t) Bn+ns (4) 
0 

where By,4n are the Bernoulli numbers. Note that the right-hand side of (4) 
is zero, when m+n is odd. Also all Bernoulli numbers 6, with indices odd 
and larger than one, are zero. 


Proof. For the proof, see [5]. 


We can express the standard basis {1,t,t?, ... ,¢’"} in terms of the Bell 
polynomials as follows: 


— 3Bo(t) + Ba(t), 
+ 11Bo(t) — 6B3(t) + Ba(t). 


e™ = S0(-1)"" [m,n] Bn(t), (5) 


n=0 


where [m,n] are called Sterling numbers of the first kind and given as follows 
(see [22]): 
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m—1%,-1 im—n—-1—1 
[m, n] _ ; 41 tot bn = ; ; ce ; OL eee eee 
OSt1<++<tm-n<m 71=0 22=0 tim—n=0 


We can write (5) as in the following way: 
t? = d;B,(t), j =0,1,2,...,m—1, 


where dj, j = 0,1,2,...,m—1, are constant vectors of order 1 x (m+ 1) 
and are given as 


dy = )(-I* Gd. 


Some properties of Sterling numbers of the first kind are as follows (see 
[4, 23]): 
(1) m, 0] = dm,03 


(2) [m,1] = (n— 1)! 
(3) [m,n] = (m—1)[m—1,n] +[m—1,n—-1], 
(4) [m,m] = 1, 


(5) [m,n] =0,n >m. 


Some sterling numbers of the first kind are given in Table 2: 


Table 2: Sterling numbers of the first kind 
[m, 1] [m,1] [m,2] [m,3] [m,4] [m,5] [m,6] [m,7] [m,8] [m,9] 


m 
0 1 

1 0 1 

2 0 1 1 

3. (0 2 3 1 

4 0 6 11 6 1 

5 (0 24 50 35 10 1 

6 0 1200 274 225 85 15 1 

7 0 720 1764 1624 735 175 21 1 

8 0 5040 13068 13132 6769 1960 322 28 1 

9 0 40320 109584 118124 67284 22249 4536 546 36 1 


2.2 The operational matrix of the Bell polynomials 


In this section, we define the operational matrix of derivative for the Bell 
polynomials. Suppose that ®,,(to) on [a, 8] is defined as 
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®,,,(to) = [Bo(to), Bi(to),..., Bm(to)]”. (6) 
Then we have 
®;,,(to) = De Pm (to), (7) 


where Dg is the (m+ 1) x (m+ 1) operational matrix of derivative for the 
Bell polynomials and it is given as follows: 


0 0 0 0 0 0 
1 0 0 0 0 0 

i eo & @ 0 

1 3 3 0 0 0 

Dg=| 1 4 6 4 O 0 
1 5 10 10 5 0 

: : 0 

GIA) (Gee Ga) 0 


Theorem 2. For all i,j €¢ Z* and i < j, we have 


Bi (to) = s (7) Bi(to). 


i=0 


Proof. Using (2) and (3), we have 


then j 
Balto) + Bi(to) = Yo (7) Bue) 
Therefore 
Bi(to) = y (7) Bstto) - By (to) = y (7) Bitte) + (2) By (to) — Bio) 
so that 
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2.3 Function approximation 


Suppose that H = L?({a,8]) and that Y = span{Bo(t), Bi(t),..., Bm(t)}. 
Since Y is a finite-dimensional and closed subspace, so Y is a complete sub- 
space of H; see [12]. Thus for an arbitrary element h(t) in H, there exists a 
unique best approximation out of Y such as ho(t), such that 


for all g(t) € Y, || A(t) — holt) IIls|| kh) — 9) |] - 


Since ho(t) € Y, there exist the unique coefficients co,...,¢m, such that 
h(t) = ho(t) = 7 cj B;(t) = CTFn(Z), (8) 
j=0 


where C' and ©,,(t) are given by 
(to) => [Bo(to), By (to), ass , Bm(to)]", C = [co, CO;-:- Poles 


and T' denotes the transposition. For computing the coefficients c;, we let 


= og BD, Ba de= etd, 01D) 2a... (9) 
: 0 
where 
ri = (Bilt), BO) = f By(Q Bitar, 
and (-,-) denotes the inner product. Let 
R= [Peale scm) K= (ki, ko, seey heals 


The matrix R is called the dual matrix of ®,,, which can be obtained as 
follows: 


B Bt ] 
Risigj+t =< Bi, B; >= f B,(t)B;(t)dt =| Sofi,p}t? So {3, qhttat 
a a p=0 q=0 


On the numerical solution of optimal control problems via ... 205 


prqtl _ qgptaqtl 


Ae aeae 
=> >t eHia} pagel 


p=0 q=0 


From (9), we have 


K’T=CTR SCT H=K TR. 


Therefore, any function h(t) € L?([a, 6]) can be expanded in terms the Bell 
polynomials as 
h(t) ~ CT ®a(t), (10) 


where 
CT = (h(t), ®,(t)) Rot. 


Also, the Bell polynomials form a complete basis over the interval [qa, 6]. 
Because using (1), we write 


®,,(t) = SX), (11) 
where 
©) =| bib bei, KOS Le aw tf E, 


and S is a lower triangular matrix with nonzero diagonal elements. Since 
from [15] the matrix S is nonsingular, hence S~! exists. Let 


D2 i Bi(t) = CT Fn (t) = 0, (12) 
j=0 

where C = [co,co,---,@m]7. using (11) and (12), we have 
CT@,,(t) = CT SX(t) =0. (13) 


According to the fundamental theorem of algebra, X(t) = {1,t,t?,...,t™} 
forms a basis; see [12]. Hence they are a linearly independent set. Therefore 
in (13), we have 


CPS =0 +C? = (eo, c9,-++,¢6m] = 0; 


hence ©®,,,(¢) is a linearly independent set. Therefore using (10), the Bell 
polynomials form a complete basis in [a, /]. 
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3 Problem statement 


In this section, we discuss the numerical solution of the OCPs. Consider the 
nonlinear system [17] 


where [A(t)],,.7, and [B(t)],,,., are matrix functions and [z(t)],,.., and 
[u(t)],,,.7 are state and control vectors, respectively. The purpose is finding 
the optimal control u(t) and the corresponding state trajectory x(t), which 
satisfy (14) and (15) while maximize (or minimize) the quadratic performance 
index 


1 


B 
Z= 527 (B)G(t)x(8) iS i (x (t)O(t)x(t) + u" (t)R(t)u(t))dt, (16) 


where [G(t)] ym and [Q(t)],..m are Symmetric positive semi-definite matri- 


ces and [R(t)],.,, is a symmetric positive definite matrix. 


4 The proposed method 


In this section, we discuss a new approach to solve the OCPs by using Bell 
polynomials. Let 


i(to) ~ @F (to) X’, 1=1,...,m, (17) 
tig(to) OF Go)? FH 1jae 5h, (18) 
where ®,,,(to) is defined in (6) and [X*] (m+t)x1 and LO ees Sac are state 


and control coefficient vectors, respectively. Then by using (7), we get 


Using (17) and (18), we have 


(to) = [OF (to) X]? = [D> Billo) XP,---) > Billo) XP]", (19) 
i=0 1i=0 
u(to) = (OF, (to)UT” = [) 2 Bilto)U},-.-. > BiltoUZ1, (20) 
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where X = (X7)(m41)xm and U = (U?)(m+i)xn are state and control coeffi- 
cient matrices, respectively. The boundary conditions in (15) can be rewritten 


as 


a(a) = 2° = d°@EG,, (to), 
x(B) = x1 = d'@EG (to), (21) 
where eee ieee and FE = [1,...,1]1x(m41) are constant vectors, and 


the symbol ® denotes the Kronecker product [13]. If x(a) or x(8) is unknown 
in (15), then we put 


x(a) = [67,(a)X]? = [> Bila)X},..., 7 Bila) XP)’, 
j i=0 


2(B) ~ [®7,(B)X]? = [)) By(B)Xj,--.. 5 BB) XP". (22) 
j=0 


4.1 Performance index approximation 


By substituting (19), (20), and (21) in (16), we obtain 


max(min)Z = 50°" G(B)x° # ail / ° BF (t)O(t)®,, (t)dt]|X 


a 


B 
+ SUTEf OL ORO m(taeu. (23) 
Let 
B B 
a / ®! (t)O(t)®,,(t)dt and P* = / ®! (t)R(t)®,,(t)dt. (24) 
By replacing (22) and (24) in (23), we obtain 
2|X,U|= aXT(P EO 4: SUTPLU, (25) 


where 
P* = 7, (8)G(8)®m(8). 


The boundary conditions in (15) can be expressed as 


qk = R(a) — xR ’ k=1,...,m, (26) 


gi =ax(b)—a} , k=1,...,m. (27) 
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We now find the extremum of (25) subject to (26) and (27) by using the 
Lagrange multiplier technique. Let 


2[X,U,d°, \1] = Z[X,U] + dA°Q° + A1Q!, (28) 


where Q° = [qj], and Q' = [qj], are constant vectors. The necessary 
condition for the extremum of (28), is 


V2Z(X,U,A A] =O: 


5 Convergence analysis 


In this section, we present the convergence of the proposed method. 
Theorem 3. Suppose that the function h : [a, 6] > R is m+ 1 times con- 


tinuously differentiable (i.e., f € C™*"[a, 8]) and that S,,=span{®,,,(t)}. If 
C’B is the best approximation of h out of S,,, then 


U Bamts — g2mt+s 
T 

< 
|h — C° Bll z2[0,6) < (m+ 1)! 2m+3 ; 


where U=max |h(™*)(t)], t € [a, 5]. 


Proof. We know that, the set {1,t,t?,...,t™} is a basis for polynomials space 
of degree m. Considering 


yi(t) = h(a) + th'(a) +-++4 —hM"(a), 


and applying Taylor expansion, we have 
t(m+1) 


a(t) — wa ()] = IG) 


|, & € (a, 6). (29) 


Since C7’B is the best approximation of h out of S,, and y;(t) € Sm, using 
(29), we obtain 


Bp 
|b — CT BI Baie.g) < h — 4120.0) = / Ia(t) — ya(t) Pde 


-{ \h! ny i7—— ie se ;)?dt 
(m+ 1)! 


B e 2m+3 _ ,2m+3 
a pra, = (Sp a 
(m+)! Ja (m+ 1)! 2n+3 
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U0 B2m+3 — q2m+3 
h—CT Bllratoa < os 
ln —C Izl0.81 S Gay am +3 


Also, we can write 


U y at ai U B2m+3(1 — (g)?m4+3) 
(m +1)! Im+3 (m+)! Im+3 - (30) 


Therefore, using (30), we have 


. B2m+s 1 — (@)2m4+3 am 3 
[kh — C7 Bl r2[0,6] < . eS a : ' 
(m+ 1)! 2m +3 (m+1)!V 2n+3 


which shows that the error vanishes as m — oo. 


sO 


[22 


6 Illustrative examples 


To illustrate the proposed technique, we solve some numerical examples and 
make a comparison with some of the results in the literature in this section. 
We implemented our method with MATLAB. 


Example 1. This example has been taken from [25]. Consider 
1 
ina | (322 (t) + u2(t))dt, (31) 
0 


subject to 
a(t) = a(t) + u(t), «(0) = 1. 
We solve the problem (31), for m = 3. Let 
Lapp(t) = XO3(t), Uapp(t) = U s(t), 
we obtain 


905 285 25 = 641 1835 195 5 
236’118’ 118” ~  * 236’ 236’ 597118" 


X =(1, 


Then, the approximate solutions are as follows: 


905 285 35 
=F Bak 
Tapp(t) = Bo — x36 Bi + 77982 — 77g Bs 
4 
= 35 13 202 nes (32) 


118 59 236 
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641 1835 195 35 
appl t) = B B B B 
Uapr(t) = — 536 Bo + 53g Bi — Gq Be + 77g P3 
35 285 1125 641 
= — +7 2 4 t 
iis’ ~ 118’ * a36'~ 236 oe 
The analytical solutions are as follows: 
Be7* ot 1 2t 
exa(t) = ———— ———e ™, 34 
Texa(t) 3e-4 41° ot get co 
Be~* oo 3 =i 
Uexalt) = 3e-4 41° 3eF+1° oe 


and the optimal value of objective function is Ze, = 2.791659975310063. 
In Table 3, we compare the absolute errors of objective function Z of 
the proposed method with methods of [25, 24] for different values of m. 
Obviously, the estimated results of our proposed scheme are coincided with 
the results in [25, 24]. Figures 1 and 2 show, respectively, plots of errors for 
state and control functions for Example 1, by using the presented method. 


Table 3: Comparing absolute errors of Z of the proposed method with the methods of 
[24, 25] for different values of m for Example 1. 


m proposed method model of [25] model of [24] 
5 1.9512 x 10-1 1.9512 x 10- 1.9512 x 10-™ 
6 1.3590 x 10719 1.3590 x 1072 1.3590 x 10710 
7 3.7582 x 107!” 3.7582 x 10712 3.7582 x 10712 
8 1.9512 x 10715 2.0000 x 107! ~—- 2.0228 x 10715 


x10 
: 


Figure 1: Absolute error for the state function «(t) for Example 1, m = 8. 
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x107 
T 


2.5 


Figure 2: Absolute error for the control function u(t) for Example 1, m = 8. 


Example 2. This example has been adapted from [19]. Consider 


1 
nine | (xo(t) + u2(t))dt, (36) 
0 
subject to 


x1 (E) => xa(t), x1(0) = 0, x1(1) = 1, &(t) = u(t), x2(0) =), 


We solve the problem (36), for m = 3. Let 


Capp(t) = 3 (E)X1, Lapp(t) = O3 (t)X2, Uapp(t) = 03 (t)U. 


we obtain 


5 1 9 3 


Veta ee ~ (3,-3.0, 0]. 
ilt (0, 599) 5h 2 [0,5 579), U [3, 3, 0, 0] 


Then, the approximate solutions are as follows: 


5 1 3 1 
1 2 3 
Se Boh e Pee 
Cpe) 5 PA + 3Bo 9°33 = 5 5 
9 3 3 
2 2 
a) => g Pi = 3 P2 = 3t— 3° 5 


Uapp(t) = 3Bp — 3By = 3 —3t. 
The analytical solutions are as follows: 


3 1 3 
l @) = ?-— =, 2? (t) =3t-— <0? =3-3t 
x ( ) 2 a Pt Tonal ) 3 9 ? Wenalt) 3 3 ’ 


212 M. R. Dadashi, A. R. Haghighi, F. Soltanian and A. Yari 
and the optimal value of objective function is Zerq = 4. 


Example 3. In this example, extracted from [10], the vibration of a spring- 
mass-damper system subjected to an external force is considered. In particu- 


x(t) 
ee 
ANNA 


ou 


F(t) fe Fit) 


(a) (b) 


Figure 3: (a) Schematic of the forced-mass-damper system assuming no friction the 
surface and (b) free body diagram of the system of part (a) for Example 2. 


lar, the response to harmonic excitations, impulses, and step forcing functions 
is examined. In many environments, rotation machinery, motors, and so on 
cause periodic motions of structures to induce vibrations into other mechan- 
ical devices and structures nearby. On summing the forces, the equation for 
the forced vibration of the system in Figure 3 becomes. It is common to 
approximate the driving forces, F(¢), 


mi(t) + ca(t) + ka(t) = FO), 


where m,c, and k are fixed numbers and F(t) represents the control force 
derived from the action of one actuator force represented by F(t) = bu(t), 
where 6 is a fixed number. The linear regulator problem specific application 
in vibration suppression, the performance index for this problem is 


zZ=5 i, " (a2(#) + au2(t))dt, (37) 


to 


where tp and ty are initial and final times, respectively. We introduce the 
usual state variable notation 


Then 


£1(t) = x(t), X(t) = X(t) = Ly 24 


The boundary conditions for this example are considered as 
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x(0) = «4(0) =1, «(1) = #(1) =0 or x1 (0) = x2(0) = 1, 21 (1) = re(1) = 0. 


We solve the problem (37), using Bell polynomials for m = 3, to = 0,t- = 
1c=2,anda=b=m=k=1. Let 


Bioslt) = Og ()Xiy Ban t) = Oy Xo, tape (t) = Oy OU. 
We obtain 
X, = {1,12,-14,3], X2 = [1,-19,9,0], U = [—7, —8,4, 3]. 


Then, the approximate solutions are as follows: 


ai a(t) = Bo + 12B, — 14By + 3B3 = 313 — 5t7 +441, 


2 o(t) = Bo — 198) + 9Bs = 94? — 10¢ +1, 


Lapp 
Uapp(t) = —7Bo — 8B, + 4B2 + 3B3 = 3t? + 13t? —t —7. 


The analytical solutions are as follows: 


1 —[—** e- isz Toa! 2 
STA agg 151 yet og ST 
r 1793 © 493 © 167°”? 
2838 _ 167, 12872 r67 67 
2 = — 1674 167 4 lg 
Pewalt) [Erg ©? — G93 6? Isin(zezt) 
| 102d age 51 ty oo OL 
rson © 423 OS T67 ”? 
2931 sr, , 15449 sory, 67 
tera(t) =[- Ta35 351 0 18in(gg7?) 
(2591 yepe 5 2826 18 cos( 67 
1263 ° rp Malls es 10 


and the optimal value of objective function is Zerg = 13.00484054100315. 
Table 4 shows the errors of x1, 72, u and approximate values of objective 

function Z, and Figures 4, 5, and 6 show the absolute errors for states of x1, #2 

and control of u functions for Example 3, by using the presented method. 


Table 4: Approximations of Z and errors of +1, #2, u, for Example 3. 


Zapp errors of x1 errors of x2 errors of u 
13.01222484276730  1.9070x 10-2  1.2649x 10-2 +1.2145x 10-1 
13.01900000000000 7.6740 x 10-® = =7.2401x 107® = 9.4197 x 1074 
13.00484749962340 2.9364x 10-8 2.6746x10-5 4.1151 x 10-4 
13.00484741498875 1.9166x 10-8 2.2206x10-® 8.6153 x 10-6 


NQoesls 
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3 x10® 


0.5 - 7 


Figure 5: Absolute errors of x2 for Example 3, m = 7. 


Example 4. This example has been taken from [2]. Consider 


an ; | (22(t) + u2(t))dt, 


subject to 
“(t) = ta(t) + u(t), (0) =1. 


In Table 5, we compare the results of the optimal values of objective func- 
tion Z of the proposed method with methods of [2, 3] for m = 8. Obviously, 
the estimated results for the minimum values of Z of the proposed method 
are more accurate than the results in models of [2, 3]. Figures 7 and 8 illus- 
trate the behavior of state variable x(t) and control variable u(t) respectively, 
by using the presented method for m = 8. 
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Figure 6: Absolute errors of u for Example 3, m = 7. 


Table 5: Comparing approximate values of Z of the proposed method with the methods 
of [2, 3] for Example 4, m = 8. 


m proposed method model of [2] model of [3] 
8 0.48426769622877 0.484268 0.48427 


0.9 


0.85 


Figure 7: Behavior of z(t) for Example 4, m = 8. 


Example 5. This example has been adapted from [26]. Consider 


jeep ; | (22 (t) + 23(t) + u2(t))dt, 
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0.24 


0.44 


u(t) 


-0.6+ 


-0.8+ 


Figure 8: Behavior of u(t) for Example 4, m = 8. 


subject to 


41 (t) = —a1(t) + vo(t) + u(t), «(0) = 1, 
&o(t) = —2x9(t), 


The analytical solutions of this problem are as follows: 


= 
w(t) = Se + 2.48165e~¥! + 0.018352e¥*, 


1 
ut) = 56° — 102792267 Y7" + 0.044305e%", 


and Z = 0.43198832823734928111. In Tables 6, 7, and 8, we compare the 
results of the absolute errors of x1(t), x2(t), and u(t) of the proposed method 
with the method of [26] for different points of time and values of m, respec- 
tively. Obviously, the estimated results of the proposed method are more 
accurate than the results in [26]. Figures 9, 10, and 11 illustrate the behavior 
of state variables x1 (t), 72(¢) and control variable u(t) by using the presented 
method for m = 8, respectively. 
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Table 6: Comparing absolute errors of #1 (t) of the proposed method with the method of 
[26] for Example 5, different values of m. 


proposed method model of [26] 

t m=A4 m=6 m=A4 m=6 

0.1 1.5215x 10-4 2.4791 x 10-8 6.2840 x 10-5 6.9470 x 10~§ 
0.2 2.9911x 10-5 4.1683 x 10-7 3.8617 x 10-4 5.6627 x 10-6 
0.3 1.2214x 1074 = 1.5184 x 1076 4.7516 x 107* ~—- 1.4208 x 1076 
0.4 1.0726 x 1074 ~—-.2.1591 x 1076 2.5623 x 1074 = 5.4834 x 107-7 
0.5 2.2672 x 10-5 8.9684 x 10-7 1.1356 x 10-4 —- 33.5375 x 10-8 
0.6 1.2331 x 107-* 8.4596 x 1078 4.0262 x 10-4 5.8476 x 1076 
0.7 8.1944 x 10-5 ~—- 9.0533 x 10-7 4.2904 x 10-4 3.9087 x 10-6 
0.8 8.3818 x 10-5 —- 1.8245 x 10-6 1.5915 x 1074 1.1543 x 107-7 
0.9 1.4779 x 1074 = 2.4645 x 107-7 2.1417 x 1074 1.8431 x 1076 


Table 7: Comparing absolute errors of x2(t) of the proposed method with the models of 
[26] for Example 5, different values of m. 


proposed method model of [26] 

t m=A4 m=6 m=A4 m=6 

0.1 1.6805 x 10-* 6.0534 x 10-7 5.9651 x 10-° ~—- 1.8022 x 10-6 
0.2 4.1122 x 1075 9.0967 x 1077 4.0570 x 10-4 ~—- 1.2818 x 10-6 
0.3 1.2262 x 107-4 ~=1.4305 x 107-7 5.1060 x 1074 ~—-:1.7540 x 1076 
0.4 1.1136 x 1074 = 8.0849 x 1077 2.8780 x 10-4 —- 2.1099 x 10-6 
0.5 2.4713 x 1075 ~—- 1.0834 x 1077 1.0392 x 10-4 6.6340 x 107-7 
0.6 1.3364 x 10-4 = 8.2323 x 10-7 4.2076 x 1074 =. 2.8671 x 1076 
0.7 9.2876 x 10-5 6.4668 x 1078 6.6499 x 10-4 —- 1.5552 x 10-6 
0.8 8.3272 x 1075 ~— 8.9367 x 107-7 1.9246 x 10-4 —-1.5283 x 1076 
0.9 1.5574 x 10-4 7.4550 x 10-7 1.9951 x 1074 =: 1.2855 x 107-8 


Table 8: Comparing absolute errors of u(t) 
[26] for Example 5, different values of m. 


proposed method model of [26] 

£ m=4 m=6 m=A4 m=6 

0.1 2.9430 x 1075 1.1055 x 1076 1.0466 x 1075 ~—- 2.5978 x 10-6 
0.2 9.1300 x 10-7 ~—- 1.2850 x 10-6 8.0579 x 10-5 2.1213 x 10-6 
0.3 3.4541 x 1075 5.5935 x 107-7 8.9374 x 10-5 8.6406 x 10-7 
0.4 2.7359 x 10-5 ~—- 9.289 x 10710 2.5833 x 10-5 5.0892 x 10-7 
0.5 4.9785 x 107-& = 1.4324 x 1078 6.7342 x 1075 ~—-1.3463 x 1076 
0.6 2.7466 x 1075 ~—-5.9741 x 1078 1.3411 x 1075 = 1.9603 x 10-8 
0.7 1.5686 x 1075 = 7.1574 x 1077 1.3507 x 10-4 ~—- 1.3192 x 107-6 
0.8 2.2792 x 1075 —:1.3485 x 1076 6.8786 x 107° ~—-1.1611 x 1076 
0.9 3.4034 x 10-5 ~—- 1.1203 x 10-6 1.3754 10-5 = 1.0032 x 10-7 


of the proposed method with the models of 
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Figure 9: Behavior of x(t) for Example 5, m = 8. 


Figure 10: Behavior of x2(t)for Example 5, m = 8. 


0.1 


Figure 11: Behavior of u(t) for Example 5, m = 8. 
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7 Conclusion 


In the presented study, we have provided a computational approach to solve 
linear constrained quadratic OCPs by employing the Bell polynomials. Sev- 
eral test problems have been studied to illustrate the applicability and effi- 
ciency of the proposed method. Results have revealed that the new technique 
can solve the OCPs effectively. Comparison of the proposed method with 
other previous methods showed that this method is accurate. 
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